Light paths of normal and phantom Einstein-Maxwell-dilaton black 

holes 



Mustapha Azreg-A'inou 
Ba§kent University, Department of Mathematics, Baglica Campus, Ankara, Turkey 

CM 
Q-i 

D ■ Abstract 

CO. 

, Null geodesies of normal and phantom Einstein-Maxwell-dilaton black holes are determined analytically 

by the Weierstrass elliptic functions. The black hole parameters other than the mass enter, with the appropri- 
ate signs, the formula for the angle of deflection to the second order in the inverse of the impact parameter 
allowing for the identification of the nature of matter (phantom or normal). Observation of the so called 
relativistic images on photon spheres also allows for such a determination. Scattering experiences may favor 
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black holes of Einstein-anti-Maxwell-dilatonic theory for their high relative discrepancy with respect to the 
Schwarzschild value. For the cases we restrict ourselves to, phantom black holes are characterized by the 
absence of many-world and two-world null geodesies. 



PACS numbers: 04.20.Gz, 04.50.Gh, 04.70.-s 



> 

CO 
<N 

in 

G\ " 1 Introduction 

o: 

Most experimental settings for testing gravitational theories are designed to evaluate trajectories of light rays. 
Accuracy in this field is a growing interest. From this point of view, the leading experimental settings are 
^ ■ aiming to achieve high accuracies beyond the known first order level and to reach a sensitivity of 1 part in 10 9 
in measuring the Eddington parameter y [lj, which is am important parameter in post-Newtonian formalism. 

On the theoretical front, workers have been striving to evaluate exactly light paths using (hyper)-elliptic 
functions mainly the Weierstrass elliptic functions denoted by p (HO. From the one hand, this has provided 
answers to some open questions, for instance, whether the cosmological constant could be a cause of the Pioneer 
anomaly [4 ], has raised the question to whether lensing could be used as a test of the Cosmic Censorship Q 
and has lead to discover new light paths, the Pascal Limacon trajectories for black holes with cosmological 
constant j6). From the other hand, the analytical solutions derived so far, IP7M231 to mention but a few, could 
be useful for any of the experimental settings aiming to test gravitational theories. Moreover, they provide 
new academic techniques for tackling motion of massive and massless particles in the geometries of various 
gravitational fields and may serve as references for testing the accuracy of numerical methods (241 and provide 
unique benchmarks for testing and improving perturbation and decomposition methods (25]. For that purpose it 
is very helpful to have relatively simple solutions. 
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In case of spherical symmetry, one of the equations governing geodesic motion reduces to 



dr \ 2 
d0/ 



P(r) 



(1.1) 



where P(r) is a polynomial function of the radial variable r, the parameters of the solution and the constants 
of motion. Depending on the dimension of the space-time, P(r) may be reduced, as described in [16], to a 
polynomial of degree 3 or 5. We are interested in the former case and we assume that d 1 - 1 b is brought to 



by coordinate transformations. Here g2,g3 depend on the parameters of the solution and the constants of motion. 
So far no special terminology has been introduced to simply notations and expressions. We introduce the fol- 
lowing terminology to describe (11.21 ) and the related polynomial and coordinates. We shall call (11.21 ) Weierstrass 
differential equation^, w(p) = 4p 3 — gift — g3 Weierstrass polynomial and (p, 0) Weierstrass coordinates. 

We bring dl.lb to d 1 -21) by a series of coordinate transformations relating r to Weierstrass radial coordinate p 
where p(r) is a non-trivial and nonlinear transformation; however, 0(0) is a linear transformation and in many 
cases = 0, where <p is the azimuthal angle. 

Most workers prefer to use the effective potential approach by which they determine all planar trajectories 
(absorbed paths (captured photons), scattering paths, trapped or confined paths, (un)stable circular paths, spi- 
ral paths approaching the circular paths from above and/or below and some other special closed curves). The 
method we shall apply is entirely based on the properties of Weierstrass differential equation and of its poly- 
nomial. We shall develop and use this method, which has been used in [9l[T9] (and partially used in [6j[8)), 
leading to a systematic approach for all problems governed by (1 1 - 2b - This will allows us to determine all types 
of trajectories. 

None of the papers mentioned above has ever dealt with light paths of normal black holes of Einstein- 
Maxwell-dilaton (EMD). One of the purposes of this paper is to address this question; the other one is to extend 
the analysis to that of light paths of phantom black holes of EMD and to draw a comparison between the systems 
of trajectories for a given ratio of charge to mass squared (a 2 = q 2 /M 2 ). 

In a phantom gravitating field theory one or more of the matter fields appear in the action with an unusual 
sign of the kinetic term , so that they are coupled repulsively to gravity. In the case of "phantom EMD" theory, 
which is also a short term for the theory, we may have a number of ways the matter fields are coupled to 
gravity: Einstein-anti-Maxwell-anti-dilaton, Einstein-Maxwell-anti-dilaton, Einstein-Maxwell-dilaton and so 
on. The presence of phantom fields continues to receive support from both collected observational data [26 ] and 
theoretical models |[27l . 

The static, spherically symmetric black hole solutions to EMD theory with phantom Maxwell and/or dilaton 
field were derived and their causal structure was analyzed, among which one finds nine classes of asymptoti- 
cally flat and two classes of non-asymptotically flat phantom black holes |[28l . In a subsequent work (29), these 
solutions have been generalized to multicenter solutions of phantom EMD. Recently, their thermodynamic prop- 

'This is justified since the differential equation satisfied by Weierstrass elliptic functions p(&) is just dl.2t upon replacing p by p. 




(1.2) 
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erties and stability were investigated too QUI . One of the remaining tasks is to investigate their null geodesies 
to see how phantom fields may affect the light paths, particularly the angles of deflection, the photon spheres 
and related lensing effects. Deviations of the angle of deflection from the Schwarzschild value are generally due 
to extensions in the theory (inclusion of Maxwell fields or and scalar ones, cosmological constant and so on), 
departure from spherical symmetry or motion of the solution itself (mostly rotation). In this paper we examine 
the case due to the inclusion of (anti)-dilatonic and/or (anti)-Maxwell fields. 

In Section |2] we consider the cosh and sinh black hole solutions of the generalized phantom EMD, which 
depend on three parameters (M,q, y). We evaluate, and discuss, the angle of deflection <5</> to the second order of 
approximation in the inverse of the impact parameter as function of the black hole three parameters. Figures, re- 
lying on exact formulas, depicting <50 for phantom and normal black holes are plotted against the Schwarzschild 
angle of deflection for different values of the parameters. The relative discrepancy is also discussed and plotted 
showing high values from some set of the parameters. 

In Section[3]we introduce the Weierstrass elliptic functions and use and develop the method based on Weier- 
strass polynomial to determine exactly all kind of null geodesies to any spherically symmetric geometry pro- 
vided the equation of (planar) motion of light rays may be brought to (1 1 -2b - Applications are considered in 
Sections |4] and [5] In the former section we consider the case y = 1 and show that the problem of determining 
the null geodesies of normal Reissner-Nordstrom black holes by the method based on Weierstrass polynomial, 
which was initiated in |fl9l , is tractable analytically and extend the analysis to phantom Reissner-Nordstrom 
black holes upon applying the results of Section [3l and in the latter one we consider the case y = and deter- 
mine all the null geodesies of the phantom cosh and normal sinh EMD black holes by mere comparison to the 
work done in Section [3j In Sections [3] to [5J we do not aim to go into the details of each null geodesic motion; 
rather, we present a general procedure (Section[3]) by which we discuss some type of null geodesic motions, the 
nature of existing divergencies and present exact reference and standard formulas for specific geodesies and for 
the angle of deflection. The paper ends with a Conclusion Section and an Appendix one. 

2 The deflection angle of light paths in the cosh-sinh solutions of EMD 

The action for EMD theory with phantom Maxwell and/or dilaton field reads 

S = -J dV^S [^-2T Jl ^%<p5 v <p + T] 2 e 2A % v F'' v ], (2.1) 

where X is the real dilaton-Maxwell coupling constant, and 7]i = ±1, T]2 = ±1. Normal EMD corresponds to 
T]2 = l)i = +1, while phantom couplings of the dilaton field cp or/and Maxwell field F = dA are obtained for 
Tji = — 1 or/and 772 = — 1. 
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The metrics of the so called cosh and sinh solutions, derived in [28], take the form 



d, 2 = f+ f_dt 2 - u'/zw- - r y_-?dtf (2.2) 

F = -*drA*, ^ = fs\ /± = 1-^, r^ 1 '^ 2 



r ' ' 1 + T7xA 2 

tj 2 (1 + t] 1 A 2 ) <0 for cosh, tj 2 (1 + tjiA 2 ) >0 for sinh (2.3) 

/€ (-»,-l) u[l,+~) if T] 1 = -l (2.4) 

ye(-l,+l]if T]! = +l (2.5) 

where we have introduced the parameter 7 following the notation of OUll 2 ! 

These are asymptotically flat spherically symmetric black holes of mass M, electric charge q and event 
horizon r + > related by |[28l 

M= r++ J r - (2.6) 



<Z = ±y^T] 2 r + r_ (2.7) 

where we have substituted, into the original formula ofq, I + 771A 2 = 2/(1 + /). Since q is real, r_ and 772(1 + 7) 
must have the same sign. Using this fact in (12.3b . we have r_ < for the cosh solution and r > for the sinh 
one. 

As shown in subsection 4.1 case 2. (d) of [28], r = corresponds to a singularity for the cosh solution where 
geodesies terminate (a Penrose diagram is given in figure 1 of [28]). Similarly, in subsection 4.3 case 1. (a) (ii) 
of |[28l it is established that, for generic values of 1 + 771 A 2 as this is the case for 7 = (to which we restrict 
ourselves in Section [5]), r = r_ is a null singularity for the sinh solution (a Penrose diagram is given in figure 3 
of |[28ll ). The curvature scalar of (12.21 ) diverges at these two points for 7 = 

2r 3 (r — r_) 2 

Expressing (r + , r_) in terms of M and a 2 = q 2 /M 2 , one obtains 

r+ = 2M , r_ = T7 2 Ma 2 , if 7 = (2.8) 



M-Ji 2rj 2 M 2 a 2 yy / 2rj 2 Ya 2 

r + =M + ^,r_ = = ,' 2 . , M = MA\- f : 2! . , Vy>-1. (2.9) 

7 (1 + 7K V (1 + 7) 



[The limit 7 ->• in leads to K2M \. 

2 For the sinh solution the case 772(1 + T71 A 2 ) < 0, which would lead to r_ < 0, is not possible 1281 
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Related to the two Killing vectors (d t ,d<j,) are the two constants of motion (E,L) given by 



r 2 sin 2 0/ i-r# =L 



(2.10) 
(2.11) 



Since (12.21 ) is endowed with spherical symmetry, the motion happens in a plane through the origin. Letting the 
plane be 8 = 7i/2 and inserting (12.10||2~TTT) into the line element (12.21) : 

(with e = 1, for massive, massless particles, respectively) we bring it to 

L 2 



r 2f_ 



For scattering states E 2 — e > 0. Eliminating z in (12.1 II 12.121 ) we arrive at 



(2.12) 



dr 
d0 



4 r z 



L 2 



E 2 ~Ufl 



which we re- write in terms of w = 1 jr as 



.2(1-7) 



J L 2 



E 2 -f + fl 



£+■ 



L 2 u 2 



f. 



1-7 



(2.13) 



From now on we take e = so that the condition for light scattering is E 2 > 0. Now, let g(u) be the function 

g(u):=L 2 u 2 f + f 2 J- 1 

and u n = \/ r n be the point on the light scattering geodesic nearest the origin where gx (w n ) = 0. Since E 2 = g(u n ), 
the angle of deflection, which is twice the variation of <p minus %, takes the form 



8<j)=2L 



Setting u = u n x, we bring (12.141 ) to the form 



f_ r ( U )y/g( Un )-g(u) 



n. 



(2.14) 



50 =2 



Lm„\/1 — x 2 



Ax 



o fl-\u)^/g{u n )-g{u) VT^x 2 " 



(2.15) 
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With u n <C 1 for scattering with large values of the impact parameter (b = L/E), we obtain 



Lu n VT 



f- Y (u)Vs( u n)-g(u) 



1 + 



[(y-l)r_+2M] 1-x 3 



1-x 2 



(7 — l)r_x 



+ 



8(1 +x) 2 



3r^.(l+x + x 2 ) 2 + 2r_r + [27-l+27x + (67- l)x z + 2x 3 +x 4 ] 

+ r 2 [4T 2 - 1 + 6(27- l)x + (8T 2 + l)x 2 + 6x 3 + 3x 4 ] J u\ + 0[w„] 



where we have used ( 12.61 ): r + + yr_ = 2M. Performing the integrations over x we obtain 



50 



AM 



+ { -2M[r + + r_(2 7 - 1)] + -^[15r 2 +6r_r+(4 7 - 1)+^ (I67 2 - 1)]} \ + 0[1 / r n ] 3 . (2.16) 

16 r„ 



Since the values of 7 depend on 17 1 according to (I2.4L 12.51 ) and the sign of r_ is that of 172(1 + 7) by (12.7I ). 
the deflection angle depends on the type of EMD under investigation. From (12.81 12.91) one sees that, for both 
cases 7 = and 7 7^ 0, the limit case q = corresponds to r_ = and r + = 2M. Thus in the limit q — > 0, 50 
approaches the value 50 (r_ = 0,r + = 2M), which is the Schwarzschild angle of deflection 50s: 



lim 50 = 50(r_ = 0,r + = 2M) = 50 s = ^ + il^l^ i + 0[ i/ rJ 3 . 

Using this along with (12.8112.91 ) in (12.161 ) we express 50 in terms of the charges (M,q) and 50s 

1 



(2.17) 



50 = 50s 



7lM 2 a 2 



16 



r 4(3tt-8) 21 
[772^— '- + a 21 



+ ■ 



if 7 = 



5 $-8^ T7 ^ (y-l)[16y-^(y+l)]M|M-^| + ^y(7y-l) g 2 1 1 | 



4M 



+ 



■2M(M + ^#) + tj 2 



87 2 

4(l-27)a 2 M 3 
(l + y)(M + 



Vy>-1 



(2.18) 
(2.19) 
(2.20) 



7T 

+ 16 



15(M + ^) 2 + 



4(16t 2 - l)a 4 M 4 12(47- l)a 2 M 2 
+ 



(1 + 7) 2 (M + ^) 2 



1 + 7 



Vy>-1 



where we have made use of M - = T7 2 |M - Ji\ and ( 12.9I ). [The limit 7 ->■ in (12.191 ) or in (12.201 ) leads 
to (12.18l )l. Thus to the first order of approximation in l/r n all normal and phantom black holes deflect light paths 

in the same way with 50 = 4M/r n H . To the second order of approximation in l/r„, the added contribution 

to the Schwarzschild one (second terms in (12.181 12.191 )) does not depend on the sign of q but depends on the 
signs of (T]i,T}2). 

First consider the special case 7 = 0, which corresponds to T71 = +1. In normal EMD (T72 = +1) we have 
50 < 50s- In phantom EMD (172 = —1), which is E-anti-MD theory, we have 50 > 50s provided we restrict 
ourselves to the physical case a 2 < 1 (4(371 — 8)/7T ~ 1.8). Thus, in the presence of phantom fields, light rays are 
more deflected than in the normal case. Phantom fields causes light rays to bent with an angle (37C — 8)g 2 / (2r 2 ) 
larger than the angle of deflection caused by normal fields. This difference is independent on the sign of q but 
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Figure 1: The relative discrepancy Rg^ = (3<j> - 8<p s )/d<j> s vs. y for fixed (M=l,a 2 = 1/16, u„ = 0.05) for E-anti-M-(anti)-D. R Sl p 
is not always negative for y < 0: it vanishes then becomes positive for some yo between —0.1 and —0.05. 
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Figure 2: The relative discrepancy = {dip - 5<ps)/8<f>s vs. y for fixed (M = I, a 2 = 1/16, «„ = 0.05) for EM-(anti)-D. Rg^, is not 
always positive for y < 0: it vanishes then becomes negative for some yo between —0.1 and —0.05. 



depends on the mass of the black hole through r n . Using (12.171 ) we obtain 

■M cr ~ — — — a z (50 s ) 



2r 2 



32 



(2.21) 



which is 0.22% of the Schwarzschild value 8(j> s if Af = 1, a 2 = 1/4, u n = 0.05 and 0.89% of it (the exact value 
is 1.03) if Af = 1, a 2 = 1, u„ = 0.05. 

The case y = 1 is phantom or normal Reissner-Nordstrom black hole. With 



_ 3KM 2 a 2 1 

50 = 50 S -T] 2 — + if 7=1 



(2.22) 



we confirm the previous conclusions: 8(j) < 8(ps for normal Reissner-Nordstrom black holes and 8(j) > 8(f)s for 
phantom ones. A phantom Reissner-Nordstrom black hole deflects light with an angle 3nq 2 /(2r 2 ) larger than 
the deflection angle caused by a normal Reissner-Nordstrom black hole 



2rl 32 V ^' 



(2.23) 



which is 1.47 % of the Schwarzschild value 8<j) s if M = 1, a 2 = 1/4, u„ = 0.05 and 5.89 % of it (the exact value 
is 6.56) if M= 1, a 2 = 1, u n = 0.05. 
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Figure 3: The relative discrepancy R^^ = (50 — 5<ps)/8(j>s vs. y for fixed (M = I, a 2 = 1/16, w„ = 0.1) for E-anti-M-(anti)-D. Rg^ is 
not always negative for y < 0: it vanishes then becomes positive for some yo between —0. 1 and —0.05. 
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Figure 4: The relative discrepancy Rg^ = (8<p - S<j> s )/5<j> s vs. y for fixed (M = I, a 2 = 1/16, u„ = 0.3) for EM-(anti)-D. R Sl p is not 
always positive for y < 0: it vanishes then becomes negative for some yo between —0.15 and —0.1. 



Now, consider the case 7^0 and 7^1. Here again we confirm the previous conclusions: 50 < 50s for 
normal black holes and 50 > dips for phantom ones provided y is large enough. This is no longer true if y is 
closer to —1 as the coefficient of \/r 2 becomes too large invalidating (12.201) . Let us look at the derivative of the 
relative discrepancy function Rg^ 

8<p - 50s 



R 



Sip 



50s 



dRst _ 2LE 2 
dy 



ln(/-)d« 



(2.24) 



(2.25) 



50s Jo /l~ r ( M )[£2_ g ( M )]3/2 ' 

which is positive for phantom black holes (/_ > 1) and negative for normal ones (/_ < 1). Thus, RsqiY) (all 
other parameters being fixed) increases for phantom black holes and decreases for normal ones. Figures (Q] 
to |4]>), which have been plotted using the exact formula ( 12.151 ), show the existence of a zero 70 beyond which 
Rs<p > (50 > 50s) for phantom black holes and 7?g^ < (50 < 50s) for normal ones. 

Figures dT] to |4]> have been plotted for fixed (M = I, a 2 = 1/16) and relatively large values of r n (u n = 
0.05, 0.1) or small values of r n up to the photon sphere (u n = 0.3). Based on these and on some other figures (not 
shown in this paper) for larger values of a 2 up to 1 and different values of u„, we can draw a general conclusion: 
All the other parameters being fixed, there is always a root 70 in the interval (—0.2,0) to Rsq(y) = 0. Said 
otherwise, for some values of 7 in the interval (—0.2,0), it seems there is always a critical value r n = r c , larger 
than the photon sphere, where Rg^ = 0. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.05 0.1 0.15 0.2 0.25 0.3 

Figure 5: The angle of deflection (Eq. l |2.15t ) in radians vs. u„ = l/r n (r„ is the point on the null geodesic nearest the origin). In 
both plots the intermediate graph is the Schwarzschild value 8<j>s. (a): Phantom EMD cosh (772 = — 1) and normal EMD sinh (772 = +1) 
black holes for y = > 70, M = 1 and a 2 = 1/4. (b): Phantom EMD cosh (772 = —1) and normal EMD sinh (772 = +1) black holes for 
y = 0, M = 1 and a 2 = 1 . 




0.05 0.1 0.15 0.2 0.25 0.3 0.05 0.1 0.15 0.2 0.25 

Figure 6: The angle of deflection (Eq. J2.15t ) in radians vs. u„ = l/r n (r„ is the point on the null geodesic nearest the origin). In 
both plots the intermediate graph is the Schwarzschild value 50s- (a): Phantom Reissner-Nordstrom (772 = — 1) and normal Reissner- 
Nordstrom (772 = +1) black holes for y = 1 > yo, M = 1 and a 2 = 1/16. (b): Phantom Reissner-Nordstrom (772 = —1) and normal 
Reissner-Nordstrom (772 = +1) black holes for y = 1, M= 1 and a 2 = 1. 

As 7 — > +00, R s ^ approaches the limit 

2 r l e Kx dx n + 8(j) s 

8(j> s Jo y/e K - (l-K + )- x 2 e K - x (l - K+x) <50s 
K± =Mw„(v / l-2T7 2 a 2 ±l). 

As Figures ((5] to [8]) reveal, the vertical spacing \8<j) — 8 <ps | depends slightly on 7, which itself depends on 
Tji, and increases with a 2 . In the extreme case (a 2 = 1), the winding number for phantom black holes with 
T]2 = — 1 and 7 > 70 (regardless the sign of Tji) diverges near u n ~ 0.3, a value for which the angle of deflection 
for phantom (rji = — 1) or normal (171 = +1) black holes with T72 = +1 is less than few radians. As we shall see 
later, this is a consequence of the fact that the photon sphere for black holes with r/2 = — 1 (black holes where 
the Maxwell field F is coupled repulsively to gravity) and 7 > 70 is larger than 3M, which is the Schwarzschild 
limit, allowing photons to orbit the hole at larger, ever-decreasing, radii. The Schwarzschild limit 3M is larger 
than the photon sphere for black holes with 772 = +1 and 7 > 7o- If 7 < 70, all that said for black holes with 
T]2 = — 1 (respectively r/2 = +1) applies to black holes with T72 = +1 (respectively T72 = — 1). 

It is useful to express 8<p in terms of the charges (M,q) and the impact parameter b = L/E. Solving the 
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Figure 7: The angle of deflection (Eq. J2.15I in radians vs. u„ = 1 /r n (r„ is the point on the null geodesic nearest the origin). In both 
plots the intermediate graph is the Schwarzschild value 50s ■ ( a ) : E-anti-M-anti-D and EM-anti-D black holes for y = 10 > yo, M = 1 
and a 2 = 1/16. (b): E-anti-M-anti-D and EM-anti-D black holes for y= 10, M = 1 and a 2 = 1. 
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Figure 8: The angle of deflection (Eq. l !2.15t ) in radians vs. u„ = i/r n (r n is the point on the null geodesic nearest the origin). In 
both plots the intermediate graph is the Schwarzschild value 8$$. (a): E-anti-MD and EMD black holes for y = — 1/2 < yo, M = 1 and 
a 2 = 1 / 16. (b): E-anti-MD and EMD black holes for y = - 1/2, M = 1 and a 2 = 1 . 



equation E 2 =g(u n ) for u n = l/r n to the second order, we obtain 

1 r+ + r(2y-l) 1 r . ^ 

Using this in (12.161 ) we derive the desired equation 

50 = ^ + ^-[I5r 2 + + 6r_r + (4 7 - 1) + ^(16^ - 1)] -j + 0[\/b] 3 . (2.26) 
b 16 o z 

The Schwarzschild value, corresponding to the limit q — > 0, is again given by 

50s = 50 (r_ = 0, r + = 2M) = ^ + ^ + 0[1/^] 3 . (2.27) 
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In terms of Sips, as given in (12.27b . b and (M,q), the expressions (12.18ll2~T9ll2.22l ) become, respectively 

^ = ^s-^[12t] 2 + ^]1 + - ) if 7 = (2.28) 

g0 = 'ft, - rj 2 [ (y " 1)[16y - ± I ± 1 ■! + ■■., V y>-1 (2.29) 

8y 2 Z? 2 

g0 = ^ s -T 72 (37r + 8) ' ?2 l + ---, if y = 1 ■ (2-30) 

3 Determination of geodesies by Weierstrass elliptic functions: general proce- 
dure 

Any differential equation of the form (11.21 ) has a unique solution in terms of the Weierstrass elliptic functions of 
the form OUH 

p = p{& + C) 

where C is generally a complex constant. 

p(&) is an even single- valued doubly-periodic function with half periods (ft), ft)') chosen in such a way that 
Im(ft)' /ou) > 0. When the Weiersttass polynomial w(p) = 4p 3 — g 2 p — g3 = 4(p — e\)(p — e 2 )(p — £3) has three 
real roots (e3,e 2 ,ei), there are three half periods (C0\, 0>2, 6O3) depending on (ft), ft)') such that 

j&(<Bfc) = e t , (A =1,2,3). (3.1) 

In order to have < e 2 < e\ we choose the three half periods (ft^ , ft^, ft^) to satisfy (2| 

ft)i = ft) > , ft) 3 = ft)' (with - ift)' > 0) , «2 = -ft) - ft)' = — ft)i - ft>3 • (3.2) 

The expression of ft^ is a consequence of e\ + e 2 + e 3 = 0. 

3.1 Three distinct real roots 

The Weierstrass polynomial w(p) will have three real roots if 

g 2 >0 and A = ^-27gf>0. (3.3) 
We parameterize the (real) roots by the angle < f] < % as follows 




(3.4) 
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The signs of < 0, e\ > and sin r) > are well defined, < e 2 < e\ and the sign of e 2 depends on gy. 



e 2 -g3<0 and e 2 = if g 3 = (g 3 = 4<?ie 2 e 3 ) • ( 3 -5) 
Motion is possible where w(p) > 0: 

<?3 < P < <?2 or P > <?i • (3.6) 
Conversely, we can reverse (13.11 ) and express the half periods ((Oi ,0)2,(03) in terms of the rootH Q. 

^=r^==r^= (3.7) 

■M y^p) ^3 vw(p) 

« s= r * =i r_* . ( 3. 8) 

•/e 3 \/w(p) \/— w(p) 

Let poo, Po be the values of p corresponding to r = +00 r = 0, respectively, and let p + , p_ correspond 

1 



to r = r + , r = r_, if there are an)Q. A singularity is denoted by p s j ng (r s i ng ), which may be any of po, p_ 



depending on the theory. In a general physical situation (poo,Po,p+,P-, • • • ,£3, £2,^1) WQ functional of the 
vector of parameters p = (charges, constants of motion) = (M,q,...,E,L,...) so that the locations of these 
points on the p-axis change with p. We shall represent (^3,^2,^1) at the same locations on the p-axis while 
(Poo,Po,p+,P-, • • •) appear on different locations depending on p. 

In order to determine all types of geodesic motion in a given geometry, one has to consider all allowed 
possible locations of {p 00 ,p a ,p + ,p_, . . .) with respect to (^3,^21^1)- Once this is done, any geodesic motion that 
can be brought to (11.21) is integrated by mere comparison with the work done in this section. We shall provide 
some examples in this section and in the next two we apply the procedure to EMD and to Reissner-Nordstrom 
phantom and normal black holes. In order to illustrate the procedure, we shall envisage only some locations of 
(Poo,po,p+,P-, • • •) with respect to (^3,^2,^1) most of which are related to EMD Reissner-Nordstrom phantom 
and normal black holes. 

3.1.1 Scattering and trapped paths 

We consider four possible situations. 



(a) £3 < Poo < £2 < £1 < Po = Psing- As light scatters from r = +00 — y r n — > r = +°o, the corresponding point 
on the p-axis moves from p m — > e 2 — > p»o, if p(r) is a decreasing function of r, or from p^ — > — > p^, if p{r) 
is an increasing function of r. We consider the former case throughout this section, which is too going to be the 
case for the next two sections. The solution to (11.21 ) is 

p(r)=p(&((j>)+C). 

3 There is a third expression for CO3 which appears with a misprinted sign in |3l ll91 . The correct expression is: o>$ 

4 In case of wormhole solutions, one introduces p a corresponding to the radius a of the throat. 
5 Some of which may be constants as in the case of Schwarzschild solution where < p„ = — 1 / 1 2 < ei and po 



-if. 



«3 dp 
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To fix C we may assume = at p„ = e%, corresponding to r n ; or assume = at p m . The former case looks 
simpler leading to <?2 = P(C) and thus, by (13.11 ), C = u>i or C = — CO2 (p is an even function!). We choose the 
latter solution so that by (13.21) C = ft)i + 0)3 and 



p(r) = ^(0(0) + w 1 +w 3 )- (3.9) 

The angle of deflection is 

S*=2p *£=-* = 2 ( r-f^-^-n^-ir^-n (3.10) 

where we have used the second formula in (13.71) . 

The other possible motion, the so called trapped or terminating bound path, is in the region between e\ and 
Po where w(p) > 0. If the path starts from the singularity po = p s i ng (r = 0) it reaches the farthest point py = e\ 
(r = rf) and then back to po. Again choosing = at the farthest point, the solution is 

p{r) = p(0(»+C), withC = Wi (orC=-G>i). (3.11) 



(b) £3 < Poo < ei < e\ < p+ < Po < p_ = +°°. If Po is a singularity, then this case is identical to (a) with a 
scattering path from p M — > <?2 — > Poo given by ( 13.91) and a trapped path between e\ and po given by (13- lib . 

If p_ is a singularity but po is not, then there is a trapped path between e\ and p_ given by (13- 11b - 

If neither po nor p_ is a singularity, the path is a many-world periodic bound orbit lfl6l in that the path, after 

crossing the inner horizon at r = r_ , emerges in another copy of the space-time then in another copy of it and 

so on. If we choose = at p = e\, then the solution will be given by (13.111) . 

(c) e2 < Po < e\ < Pa, < p+ < p_ = +00. Since w(p) < for p G (<?2,<?i), there are no paths that can reach or 
emanate from the origin. 

There is a path that extends from spatial infinity (poo) to the inner horizon (P-). This is not a spiral path 
since the integral 

/ dp/y^p) 

J constant > 

converges (0 remains finite). If p_ is not a singularity, the path is called a two-world scattering orbit in that the 
path emerges, after crossing the inner horizon at r = r_ , in another copy of the space-time and back to spatial 
infinity. If p_ is a singularity, we have an absorbed path from spatial infinity to the singularity. The solution 
is again p(r) = p(®(<p) +C). Since there is no farthest or nearest point on the path, we choose = at poo 
leading to p„ = p{C). Using the inverse function to p, C = dp / y / w(p) and 

p(r) = p(0(0) + C), with C= / dp/7^p~). (3.12) 

(d) Po < £3 < Poo < ei < e\ < p_ = +00. Here again no paths that can reach or emanate from the origin. 

We have a scattering path from — > — > Poo and the solution is again given by ( 13.91 ) but with different p. 
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There is another path from r = r\ (p = e\) to r = r_ (p = p_). Since in this case n is finite (p„ < ei 
r! < +00), the path is a many-world periodic bound orbit if p_ is not a singularity or a trapped path if p_ is a 
singularity. If we choose © = at p = e\, then the solution will be given by (13- 11b - 



3.1.2 Absorbed and circular paths 

Absorbed paths extend from spatial infinity (poo) to the (nearest) singularity (p s i ng ). Such paths exist if both 
points poo, p s i n g are in [^3,^2] or in [e\ , +°°). [There are no such paths in the Schwarzschild case when w(p) has 
three distinct real roots.] It is clear that there are no circular paths when w(p) has three distinct real roots. 

3.2 Two distinct real roots 

The Weierstrass polynomial w(p) will have two real roots if 



This happens when one of the local extreme values of w(p) is zero. The second condition in (13.131 ) splits into 
two cases. 

3.2.1 Stable circular and bound paths: §3 = (g2/3)\/g2/3 > 

The local maximum value of w(p), which is at p max = —(1/2) ^^2/3, is zero. We have 



Since at p max = — (1 /2)y g2/3, w{p) has a local maximum, the polynomial P{r) in (11.11 ) has a local maximum 
too at the corresponding point r max . But since P{r) <x E 2 — V{r) (compare with (I2.12|[2TT3T ». the potential V(r) 
has there a local minimum. Thus 



is a stable circular path. 

Paths in the region p > e\ are periodic: They include the periodic bound and the so called terminating bound 
(trapped) and many-world periodic bound orbits |[T6ll . No matter the location of p^ with respect to e\ is, the 
equation of motion can be integrated directly. Let t = yjp — e\ and 



g 2 > and A = gl-27gl =0. 



(3.13) 




(3.14) 



P — Pmax 



ea = -(1/2V&73 



(3.15) 



k 2 =e l -e 3 = (3/2) y / g~j3>0 



(3.16) 



then (11.21 ) reads 



d@ = 



dp dt 



2{p-e^)y/p~^Ti t 2 + k 2 
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leading to 



V^ = *tan[*(0-C)] or P->- 1+ ^2-C)y ^ 



For the Schwarzschild black hole p =M/(2r) - 1/12 021, = 0, g 2 = 1/12, g 3 = 1/216, e 3 = e 2 = -1/12, 
e\ = 1/6, k = 1/2 and the second formula in (13.171 ) is just Eq. (10) of 1 19]. 

For the phantom (7] 2 = — 1) and normal (t]2 = +1) Reissner-Nordstrdm black holes we derive in Subsec- 
tion [43] from the first formula in (13.171) the following orbit 



tan 



(b — Ci lr + — r 

+ (3.18) 



which is a trapped path for the phantom black hole (r + > r > 0) and a many-world periodic bound path for the 
normal black hole (r + > r > r_). Using the double-angle formula for tan we re- write it as 

M-r 

tan(0 - C) = = , (3. 19) 

Y / 2Mr - rj2<r _ r 

which is the correct form of the misprinted Eq. (64) of lfT9li . 



3.2.2 Unstable circular and spiral paths: gz = — (g2/3)\/g2/3 < 



The local minimum value of w(p), which is at p m j n = +(1 /2) y/ gzf$, is zero. We have 

ei = ei = + \^' €3 = ~\f^~ = 71 in TO • (3.20) 



Since at p m ; n = +(l/2)y / g2/3, w(p) has a local minimum, the potential V(r) has there a local maximum 
(compare with ( |2TT2ll2TT3l ). Thus 

P=p m n=ey =+(1/2)^73 (3.21) 

is an unstable circular path. 

Paths in the regions ej, < p < e 2 and p > e\ depend on the location of p^. Here we consider the case 
e 3 < Poo < e 2 . There are two spiral paths approaching the circle p = e\ = +(l/2)yg2/3 from 1) p m or from 
2) Po (Po > ^i hi this case!). The paths end orbiting the center at an ever 1) decreasing or 2) increasing radii r 
without, however, reaching the unstable circular path at r = r\ corresponding to p = e\. The equation of motion 
can be integrated directly. Let s = \Jp — e-s and k 2 = e\—e^ = (3/2) y / g 2 /3 > as in (13.161) . Then (11.21 ) reads 

dp ds 



d® 



2\p-ei\^p-e 3 \s 2 -k 2 \ 



and &, as well as the angle of deflection, diverge as ln|p — e\\ as p — > e\, which is a general behavior in the 
strong field limit valid for all spherically symmetric solutions li 3~Tll3"2"ll . Integration leads to 

1) a/p 3 ^ = -kco\h[k(® - C)] and 2) ^/p~^l = -A:tanh[/t(@ - C)} (3.22) 
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or 

= - ijl'mlt or <--»»>. +-»*»■ 033, 

For the Schwarzschild black hole p = M/(2r) - 1/12 (15), = 0, g 2 = 1/12, g 3 = -1/216, e 3 = -1/6, 
ei = ei = 1/12 and = 1/2 we obtain the solutions (1 1) of lfl9l . 

3.3 One real root 

The Weierstrass polynomial w(p) will have one real root with multiplicity 1 if 

A = gl-T7gl<0. (3.24) 



The sign of the real root e r 



^3 [(9*3 + V3V^A) l/3 + (9g 3 - V3V^aY /3 } 



(3.25) 



is related to that of g 3 by 

e r -g 3 >0 and e r = if g 3 = 0. (3.26) 

Motion is possible for p > e r . 

Absorbed paths exist if the range e r < p < °° includes r s i ng < r < oo. In that case, we will have e r <p^< p s i ng 
and the solution will be given by (13.121 ) where the upper limit of integration "oo" is replaced by "p s i ng "- 

If Poo < e r < p s i n g, the solution is a trapped path given by (13- 11b - 

One can envisage other situations as p s j ng < e r and so on. However, we are giving examples that are more 
or less related to EMD and Reissner-Nordstrom black holes. 

The case g2 = g 3 = implies e r = (this is the only case where the three real roots of w(p) =0 are 
equal). This is no difference from the two cases discussed above for the generic case. However, if p(r) were an 
increasing function of r and p„o > e r = 0, the angle would diverge as 

re,-=0 dp' 1 

0-C= / lim 



for paths approaching e r = from the right. This is not a logarithmic behavior as the one we have seen earlier. 
These spiral paths would approach, without reaching it, the unstable circular path at p = e r = 0. 



4 Phantom and Normal Reissner-Nordstrom black holes 

It is difficult or impossible to reduce (12.131) to (11.21 ) for any y. For y = ±2 and probably for some other values, 
the problem can be tackled semi-analytically in a similar way to what is done in fll6lfT6l : Limits to the analytical 
treatment are 1) lack of "compact" solutions to the polynomial equation (or its polynomial reduced form) P(r) = 
and/or 2) lack of solutions to the generally non-polynomial equation A = g\ — 21 g\ = 0. All that said does not 
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apply to the cases 7 = 1 and 7 = 0, which can be entirely analytically solved. In this section we investigate the 
former case and in the next one we tackle the latter case. 
Setting 7 = 1 in (l2T6l [2771 [279ll27T3l) with e = we obtain 

M = r+ - r ~ , q = ±VT] 2 r_r+, r_ = lfe|r_|, r_r + = T} 2 q 2 (4.1) 



r±=M± a/M 2 - tj 2 <7 2 = M( 1 ± V 1 - T 72« 2 ) (4-2) 
, \ 2 

Ji) =f-H 2 + 2Mu 3 -n 2 flV (4.3) 



where £ = E 2 /L 2 = l/b 2 > 0. 

Let u r = l/r r be any real root of the polynomial in the r.h.s of (14.31 ): 

£-u 2 + 2Mu 3 r -ri2q 2 u 4 r =0 [■£ = r} 2 q 2 u 2 r {u r - ii-){u r - u + ) > 0] . (4.4) 

Since £ > 0, a/Z the real roots of the polynomial are either greater than w_ > u + > or smaller that m + for 
T]2 = + 1 and they are all between m_ < and u + > for r\ 2 = — 1 (Figure |9]). 

If •£ 7^ {£ > 0), we have necessarily u r / w + , u r 7^ m_ and w r 7^ 0: It is not possible to find solutions with 
£ / and r r = r+ [QU Eqs. (65, 66)]. 

If « r is a root with multiplicity 1, following the general procedure, we introduce the radial coordinate y = 
u — u r followed by z = l/y and finally 

3Ciz + C 2 Ci C 2 Ci r,-r C 2 
12 4(m — u r ) 12 4 r — r r 12 

with 

Ci = 2i< r (3MM r -l-2T] 2 ^ 2 ^) = 2u r (l — Mu r ) —4£/u r (4.6) 

C 2 = 6Mw r - 1 - 67l 2 q 2 u 2 = 5- 6Mu r - 6£/u 2 r . (4.7) 

Eqs. (1431) reduce (1431) to (PT31 ) with 

g 2 = ^-42 2 , (G 2 = W) (4-8) 
l-54£M 2 + 36^2 2 

S3 = (49) 

A = ^[M 2 (l + 36£Q 2 ) -21£M A -Q 2 {\ +4£Q 2 ) 2 }/16 (4.10) 

= -£[\6Q 6 £ 2 + (27M 4 - 36M 2 <2 2 + SQ 4 )£ + Q 2 -M 2 }/16 

= -Q 6 £{£ - £-){£ - £+) = -ri 2 q 6 £(£ - £J){£ - £ + ) 

where we have used (14.41 ) to eliminate u r from the expressions of g 2 , #3. The new parameters (£-,£ + ) are defined 
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32a 4 

In the physical case a 2 = q 2 /M 2 < 1, < £+ < I- for the phantom black hole and £- < < £ + for the normal 
one. 

The transformation (14.51 ) "splits" the point r = r r into r+ and r~ (corresponding to u~ and respectively). 
If C\ < then r+ and r~ are sent to p = p r = +oo and p = — °°, respectively, and if C\ > the latter limits are 
reversed. As we shall see in the Appendix, it is always possible to choose the real root u r so that C\ < 0: choose 
u r to be the smallest root or the largest one for the phantom or normal black hole solution, respectively. The 
points poo, po and p± (corresponding to r = +oo, r = and r = r±) on the p-axis are given by 

i 1 C 2 5 I Mu r 

P ~ = 2^-12' P ° = 12 = 12-2^-— (412) 

P± = T7 s +t| (4-13) 

4[u± — u r ) 12 

which depend on « r whose analytic expression in terms of (M,q 2 ,£) is sizable. 

4.1 Three distinct real roots for w(p) = 

For the phantom case, we derive in the Appendix (Eq. (1A.2I )) the following order relations for the p -parameters 

e 3 < Poo < e 2 < e\ < p+ < p < p_ < p,- = . (4. 14) 

The only possible paths are scattering ones from p»o to <?2 given by ( 13.91 14.51 ) or trapped ones from any point 
e\ < p < p + to the singularity at po [from any point r + <r <r\ = \/u\, where u\ < u + is the largest root of 
2Mu 3 + q 2 u 4 = (Figure |9] (a)), to r = 0]. If the trapped path starts from p = e\ its equation is given 
by (EH ED). 

For the normal case, Eq. (IA.6I ) of the Appendix reads 

Po < e 3 < poo < e 2 < e x < p + < p_ < p,- = +oo . (4. 15) 

The only possible paths are scattering ones from p»o to e 2 given by d3.9ll4.5l ) or many- world ones from any point 
e\ < p < p + to p = +oo [from any point r + <r<r\ =\/u\ to r r = \/u r where u r > u- is the largest root of 
t — u 2 + 2Mu 3 — q 2 u A = and u\ < u + is the second largest root (Figure |9](b))]. If the many-world path starts 
from p = e\ its equation is given by (13. 11114.51) . 

4.2 Two distinct real roots for w(p) — 

This corresponds to (Eq. (13.131 )) 

-^-iriiq 2 > and (i = 0, £ = I- or £ = £ + ) . (4.16) 
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4.2.1 Case £ = 0. 



As we do and explain in the Appendix we choose u r = U-, which is the smallest (respectively largest) root for 
the phantom (respectively normal) solution when £ = 0. In this case Ci < and p (r) is a decreasing function of 
r [as r — > r+ (from the right), p — > +°°]. The order relations as given in (IA.3UA.7l) read 

T] 2 = -1 : e 3 =e 2 = p 00 = -^r < ei = p+ = \ < p < p_ = +°° (4.17) 
T?2 = + 1 : Po<e 3 = e 2 = p» = -Y2 < e i=P+=6<P- = + 00 - ( 4 - 18 ) 

Since po is a singularity for the phantom black hole, there is a trapped path for this hole from p + — > po given 
by (I3TT81I3TT91) with T] 2 = -1. 

There is a many-world periodic path for the normal black hole from p + — > p_ given by (13.171 14.51) . Us- 
ing ( 14.51 ) with u = \ / r, u r = u— , Ci = — (r + — r_)/r 2 and C 2 = (2r_ — 3r + )/r_ we have 

4(p- ei ) = ^^= - M (r+ ~f 2 (T]2 = +l) (4-19) 
r — r_ 2Mr — 772*7 — r z 

then using the first equation (13.171 ) with k = ^Je\ — es = 1/2 leads to (l3.18l[3~T9l . 

Had we chosen u r = u + , instead of u r = ii-, we would reach the same conclusions concerning the nature of 
the paths. In this case, we would have C\ = 2u + {\ —Mu + ) > 0, p(r) is an increasing function of r [as r — > 
(from the left), p — > +00] and 

T] 2 = -1 : e 3 =e 2 = poo = -— < e Y = p_ = \ < p < p+ = +°° (4.20) 
T] 2 = +l: e3=^2 = Poo = --^<po<ei=p_ = ^<p+ = +oo. (4.21) 

But instead of (14.191 ) we would obtain 

r-r- 2Mr - r\ 2 q 2 -r 2 



4(p — e\) = = - t -2 and tan 



r — r _ 



r + — r (r + —r) A L 2 J V r + — r 

4.2.2 Cases £ = £±. 

Now we consider the cases £ = £±. Note that in this case g 2 , g3 are, by (14.111 ). functions of (a 2 ,r\ 2 ) only 

1 



(4.22) 



gi = ^_-T]2{q 2 £± 

l-(54/a 2 )(^ ± ) + T] 2 36(^ ± ) 

#3 = 



216 

and that £± are functions of (M,a 2 , r\ 2 ) 



-27 + 36T] 2 a 2 - 8a 4 ± (9 - 8T] 2 a 2 ) 3 / 2 
7)2 32M 2 a 6 
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leading 



Urn = +00 (t} 2 = -1) and lim £ + = (tj 2 = ±1) . 



With b = y/TJI, the last two limits are the Schwarzschild limit for the impact parameter (3\/3M) allowing 
photons to orbit endlessly the hole around the photon sphere without reaching it. 

£ = £ + . If I = £ + and T72 = — 1, gi > and g 3 < 0. A similar case has been treated in Eqs. (13.201) to (13.23I ). 
There is a root with multiplicity 2 at p = e\ = e 2 = (1 /2)y 'gift. The corresponding root m = wi is such that 
the r.h.s of (14.31 ) reads: £ + — u 2 +2Mu 3 + q 2 u 4 = q 2 (u — u\) 2 (u — u r )(u — u 3 ) where u = u 3 corresponds to 
p = 63 = —2ei. This is the case where the point M\ is on the w-axis (Figure|9](a)). The order relations are given 
in (IA.4I ): 

e 3 < poo < e\ = e 2 < p+ < Po < P- < Pr = +°° • (4.23) 
There is an unstable circular path at 



(9 + 8a 2 )(9 + 4a 2 -3\/9 + 8a 2 ) 



24V2a 



2 



corresponding to r = r\ = \/u\ (the photon sphere) wit 

V9 + 8a 2 + 3 

H = 2 M > 3M > r+ (4.24) 

and spiral paths from r = +00 (p = p^) to r = r\ (p = ej) and from r = r + (p = p + ) to r = ri (p = <?i) given 
by (13.231 ) and (14.51 ) with & = \fie[. There is also a trapped path from r = r\ (p = ei) to r = (p = po) given 
by ( 13.231) with the + sign. In the limit a 2 — ^ 0, r\ — > 3M which is the Schwarzschild limit. 

If £ = £ + and T72 = +1, gi > and #3 < 0. There is a root with multiplicity 2 at p = e\ = e 2 = (l/2)y / < g2/3. 
The corresponding root u = mi is such that the r.h.s of (14.31 ) reads: ^ + — u 2 + 2Mm 3 — g 2 w 4 = q 2 {u — u\ ) 2 {u — 
u r ){u — ut,) where u = u 3 corresponds to p = = —2e\. This is the case where the point M3 is on the M-axis 
(Figure|9](b)). The order relations are given in (1A.8I ): 

po < e 3 < poo < e x = e 2 < p+ < p_ < p r = +°° . (4.25) 

There is an unstable circular path at 



^ 2 -9)(4a 2 -9 + 3V9-8a 2 ) 
9 ~ 61 ~ 24V2V 



u\ is the largest root of £ + — u 2 + 2Mu 3 + q 2 u 4 = when Mi is on the ;<-axis (Figure|9](a)). 
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corresponding to r = r\ = l/ui (the photon sphere) wifhQ 



V9 - 8a 2 + 3 

r+<n = - M<3M (4.26) 

and spiral paths from r = +oo (p = p^) to r = r\ (p = ei) and from r = r + (p = p + ) to r = n (p = e\) given 
by ( |3.23l > and ( 14.51 ) with & = y/3e[. There is also a many-world periodic bound path from r > r + through r = r_ 
to r = r r > which emerges in another copy of the space-time after crossing r = r_ . This is also given (13.231 ) 
with the + sign. In the limit a 2 — > 0, r\ — > 3M which is the Schwarzschild limit. 

1 = 1— We have necessarily T]2 = — 1 since l_ < for normal Reissner-Nordstrom black hole. In this case, 
the r.h.s of |0]>, i- - u 2 + 2Mu 3 + q 2 u A , has only one real rooo u\ = -[V9 + 8a 2 + 3]/(4a 2 M) < with 
multiplicity 2 and two complex rooti§. Thus, the r.h.s of (14.31) is always positive with only absorbed paths from 
spatial infinity to the singularity at r = given by (13.121 ) and (14.5b - 

4.3 One real root for w(p) — 

This corresponds to (Eq. (13.241 )) 

" + 36 " 2 + 32MV + 8 ° 2>3/2 < e < 2? + ^ + 32«'/ + 81,2)3/2 ^ 

1 > 36 " 2 - 27 ~Lt 81,2)3/2 ° 

For the phantom solution (T]2 = — 1), this is the case where the point M\ is above the w-axis and M 3 is below 
it (Figure [9] (a)). The two real roots (u r < 113) of t — u 2 + 2Mu 3 + g 2 M 4 = are negative and {u\ , u-i) are now 
complex roots, so (^1,^2) no longer exist . Eqs. dA.UIA.2l ) become 

m_ < u r < Ut, < < u + 

e 3 < Poo < p+ < po < p < p,- = +°° (4.29) 



with only absorbed paths from spatial infinity to the singularity at r = given by (13.121 ) and (14.5I ). 

For the normal solution (t]2 = +1), this is the case where the point M 3 is above the w-axis (Figure |9](b)). 
The two real roots of I — u 2 + 2Mu 3 — q 2 u 4 = satisfy m 3 < and u r > m_ > is the largest one. (u\ , M2) are 
now complex roots, so (^1,^2) no longer exist. Eqs. dA.5IIA.6l ) become 

w 3 < < u+ < U- < u r 

Po < e 3 < poo < p+ < p- < p,- = +00 . (4.30) 



7 «1 is smallest positive root of £ + — u 2 + 2Mm 3 — q 2 u^ = when the point M3 is on the ;<-axis (Figure |9](b)). 
8 «l is the only real root of £_ — m 2 + 2Mm 3 + ^ 2 m 4 = when the point M3 is on the M-axis (Figure[9](a)). 

9 In this case to reduction of j43\ does not lead to jl.2j ; rather, it leads to a similar equation with an irreducible quadratic form on 
the r.h.s., a polynomial of degree 2 with complex roots. 
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The only existing paths are two-world scattering paths from spatial infinity to r = r r = 1 ju r > given by (13.121) 
and d43T ). 



5 Null geodesies of phantom and normal EMD 

In section we restrict ourselves to the case 7 = which corresponds to T]i = +1, then ( 12.31 ) implies rj 2 = — 1 for 
the cosh solution and 772 = +1 for the sinh one. Thus we will be considering E-anti-MD for the cosh solution 
and normal EMD for the sinh one. 

Instead of u = (l/r_) — (/_/r_), we use /_ as a radial coordinate. This way we reduce (12.131) for light 
paths (e = 0) to 



dtp 

where, using (12.81) . 



d ' M = [a/!-(3a + l)/! + (3a + 0+2)/_-(a + l)]/_ (5.1) 



r + 2M 2 2 „ r 2 E 2 q 4 E 2 q 4 

r_ g 2 a 2 L 2 M 2 L 2 M 2 b 2 

In the physical case a 2 > 1, to which we restrict ourselves, a is constrained by 

a > 2 if 772 = -1 , a < -2 if tj 2 = +1 (5.3) 

for the phantom cosh and normal sinh solutions, respectively. 

The next step is to introduce the variable R = 1/ (/_ — / ) where f is a zero of the fourth order polynomial 
in /_ on the r.h.s of (15.11) . We choose fo = 0, leading to 



fdR\ 



2 



) =a-(3a + \)R+(3a + l5+2)R 2 -(a + \)R 3 . 



Vd0 J 

The final steps consist in eliminating the term in R 2 and rescaling by introducing the Weierstrass coordinates 
(p , &) defined by 

4 ! /3 3a + j8 + 2 

R =-J^W P+ ^TW (5 - 4) 

dQ = -T7 2 ^ J 1/3 ; d0, (d0-d©>O) (5-5) 

and d0 • d@ > by (15.31) . The reduced equation is (11.21 ): (dp /d®) 2 = 4p 3 — g 2 P — g3 with 

4^3 1 + 2(2 + 3a)j3+j8 2 

g2 = -3 — ( } 

_ 2 - 3(5 + 12a + 9a 2 )j8 - 6(2 + 3a)j8 2 - 2j3 3 
83 ~ 27(1 + a) 2 ' ( } 
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Note that, if a is restricted by (I5.3I ). p (r) is a decreasing function of r for all 772 ■ p (r) and its inverse function 
are given by 

_ (j3 - l)r- (3a + j8 + 2)r_ _ r_ [3a + ft + 2 - 3 -^(l + a) 2 / 3 p] 

P " 3-4V3(i + a)2/3( r _ r _) and r - j3-i-3-4 1 /3(i + a)2/3 p (5 " 8) 

so that, using r_ = 772 1 r_ | and a + 1 = — T?2|a + 1|, we arrive at dp/dr = — |r_||a + 1 1 1 / 3 / [4 1 / 3 (r — r_) 2 ]. 

In the limit r — >• r_, p — >■ — 3r_(a + l)/(r — r_) = 3|r_| |a + l|/(r — r_) for all 772- Thus the transforma- 
tion (15.81) "splits" the point r_ into rZ and rZ and sends the point rZ to p = —00 and the point rZ to p = p_ = +°°. 
The points poo, po and p + (corresponding to r = +00, r = and r = r + ) on the p-axis are given by 

JS-1 ^ 3a + /3 + 2 _ j3+2 



r 3-4V3(a + i)2/3'^ 3-4V3(« + l)2/3'^ 3 -41/3(1 + 0)2/3 

and po > for phantom black holes. If (ei,e2,ej) are real, the order relations of these roots with respect to 
(p<x,,Po,p+) depend on (a,j8) = p. This will be done for each case (phantom or normal) separately. 

Ordering (p«, po, p+) is also done separately as follows. For the phantom cosh black hole we have rZ < < 
r + < +00 which leads to (p (r) is always decreasing) 

Poo < p+ < po < P- = +°° ■ (5. 10) 

For the normal sinh black hole we have < rZ < rZ < r + < +00. But since p(r) is always decreasing, if one 
moves on the r-axis along the path r = +00 — > r + — > rZ — > rZ — > 0, the corresponding point on the p-axis moves 
along the path — > p + — > p_ = +00 — y (in a circular rotation) —00 — > p . Thus 

Po < Poo < p+ < p_ = +00. (5.11) 

Let (j6i ,p2,p3, AO be the following a-functions 

j3i, 2 = -(2 + 3a) + v/3(l + a)(l + 3a) (l->-2->+) (5.12) 
1 - 18a -27a 2 ± (l+9a)-v/(l + a)(l + 9a) 



P3 < 4 8a 
in terms of which we have 



(3-»+,4->-) (5.13) 



4 V 3 Q3-j3i)Q8-fe) ..... 
(l + «)4/3 (5 " 14) 

A _ 3 97 2 _ J8 [4 + ( 1 ~ 18a - 27a 2 )/3 - 4a/3 2 ] _ -4a/3 (/? - ft)Q3 - AO 

A = §2 — z '6 3 ~~ 



(1 + a) 2 (1 + a) 

5.1 The phantom cosh black hole: a > 2, Tfe = — 1 

In this case g 2 > for all j8 > (Eq. (lOD . At < and j8 3 > 0. 
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5.1.1 Three distinct real roots for w(p) = 

This corresponds to (Eq. (13.31 )) 

0< j 8< j 8 3 (5.15) 

which leads, using (l3T4l [5T9l [57TOT > . to 

<?3 < Pec < e 2 < e\ < p + < p < p- = +°° (p»o<0). (5.16) 

To order (poo,po,p+) with respect to {e\,e 2 ,e 3 ) we may use different methods as plotting the surfaces p + — e\ 
and so on or simply evaluate the Weierstrass polynomial and its derivatives w' = 12p 2 — g 2 and w' = 24p at 
(P~,Po,p+)- For instance, w(p + ) > 0, w'(p + ) > and w"(p + ) > 0. 

This case has been treated in subsubsection 13.1. Tl case (b). Since po is a singularity for the cosh black hole, 
there is a trapped path from e\ to po given by (13.llll531l5.8l ). The scattering path from poo — > e 2 — > Poo is given 
by (HUES EH). 



5.1.2 Two distinct real roots for w(p) = 

This corresponds to (Eq. (13.131 )) 



j8 = or /3 = /3 3 • (5.17) 



For p = 0, p+ = e\ and #3 > so that #3 = (g 2 /3) yg^/3. The relations (15.161 ) are still valid, in the limit 
we have e 3 = p^ = e 2 - This case has been treated in Eqs. ( 13.141 ) to (13.171 ). Since p is a singularity for the cosh 
black hole, there is a trapped or terminating bound path from p + = e\ (r = r + ) to po (r = 0) given by (13.1711531 
with e x = -2e 3 = -2e 2 = -2p M = 2 1 /3/[3(l + a) 2 / 3 ]: 



2 2/3( 1 + «)2/3 p = 2 +tan 2 



0-C 



2 1 /3(i + (Z )i/3 



Substituting j3 = into (15.61) . then into (13.151) and the second equation (15.81) we obtain the radius of the stable 
circular path at r = 00, as in the Schwarzschild case. 

For j6 = fa, g3 < so that #3 = —(gi/tyyfgi/? 1 - This case has been treated in Eqs. (13.201) to (13.23I ). The 
relations (15.161 ) remain valid with 2ei = 2e 2 = —e 3 = \J g 2 (fi?,) /3 where the value of g 2 at j3 = /J3 given by 

ta \ - (l+9«)[l+9« 2 + VA + a(2-3y^)] 
8l[ fo>- 48-2 1 /3 a 2( 1 + a )i/3 

There are spiral paths given by (13.231 15.51 15.81 ) which approach the unstable circular path at p = e\ from 
above(p oo )/below(p ). Using the above expression for g 2 {fi-i) and substituting yS = ^83 into (13.211 ) and the second 
equation (15.81 ) we obtain the radius r\ of the unstable circular path where A = (l + a)(l + 9a) and r + = 2M 

n = 8M ' + 3a) , («>2). (5.18) 



l + y/A-a(l0 + 21a-9y/A) + y/2AJl + y/A + a(2 + 9a-3y/A) 
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The limit q 2 — > corresponds to a — > +°°. The radius r\, as given by (15.181) . decreases from (5 + V57)r + /8 to 
the Schwarzschild limit 3r + /2 as a increases from 2 — > +00. 

5.1.3 One real root for w(p) = 

This case corresponds to (Eq. (13.241) ) 

j3>/3 3 (5.19) 

leading to 

Poo < e r < p+ < po < p- = +00 (5.20) 

where e r is the real given by (13.25I ). For a > 2 it is not possible to have g 2 = and #3 = 0, so there is no solution 
e r = with multiplicity 3. Since po is a singularity for the cosh black hole, there is a trapped path from e r to po 
for the generic case j8 > /3 3 given by (I3.llll5.5ll5.8l ). 

5.2 The normal sinh black hole: a < —2, Tfc — +1 

In this case < j3 4 < p\ < p\ < j3 3 . Thus, the condition A > (Eq. (15.141 )) ensures g 2 > 0. 

5.2.1 Three distinct real roots for w(p) = 

This corresponds to (Eq. (13.31 )) 

< 15 < p 4 or /3 > J3 3 (5.21) 

which leads, using (13T41 [5T9ll57TTb . to 

Po < e 3 < poo < e 2 < e x < p+ < p_ = +00 if < j3 < /3 4 (po» < 0) (5.22) 
e 3 < e 2 < po < ei < poo < p+ < p- = +00 if j3 > j8 3 . (5.23) 

For the case < j8 < pV which has been treated in subsubsection 13 . 1 . fl case (d), the solution for the scattering 
path from p»o — > e 2 — > p»o (r = 00 — > r 2 — > r = °°) is given by ( 13.91 15-51 15 - 8b - There is another path from r = r\ 
(p = e{) to r = r_ (p = p_), which is a trapped path for p_ is a null singularity for the sinh black hole. If we 
choose © = at p = e\, then the solution is given by (I3.llll5.5ll5.8l ). 

The case j8 > j8 3 has been treated in subsubsection 13.1.11 case (c). Since p_ is a singularity, we have an 
absorbed path from spatial infinity to the singularity. The solution is given by (13.1211531 15 - 8b - 

5.2.2 Two distinct real roots for w(p) = 

This corresponds to (Eq. (13.131 )) 

j3=0, p = fa or /3=j3 3 . (5.24) 

The discussion of the case j8 = for the normal sinh black hole is similar to that for the phantom cosh one. 
What was said in the first paragraph following ( 15.171 ) applies to this case if we replace "po" by "p_", "r = 0" 
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by "r = r_" and "cosh" by "sinh". Thus, there a trapped path from p + to the singularity p_ given by (13.1711531 



For /3 = p4, gi < so that gi = — (g2/3) y / g2/3. This case corresponds to the case j8 = ($3 of the phantom 
cosh black hole; the discussion in the second paragraph following (15.171) applies and the radius of the unstable 
circular path is obtained from (15.181 ) upon changing yA to — a/A [this is also true for £2(^4) of this case which 
is derived from giifh) of the phantom cosh black hole by VX — > — VA] 

n = 8r +" +3a > («<-2). (5.25) 



1 - VI - a ( 10 + 27a + 9 y/A) - V2A J 1 - VX + a (2 + 9a + 3 VA) 

The limit q 2 — > corresponds to a — >• — °°. The radius ri, as given by (15.251) . increases from (7 + \/l7)r + /8 > r + 
to the Schwarzschild limit 3r + /2 as |a| increases from 2 — > +00 (« decreases from —2 — > —00). 

For (5 = fa, gi > so that #3 = +(g2/3) ^fgifl. In this case £3=^2 = — 2ei (Eq. (13.141 )) and all remaining 
inequalities in (15.231) are still valid. Since p is a singularity, we have an absorbed path from spatial infinity 
to the singularity. The solution is given by (13.12[|5"31l5.8l ). The value rj = r2 corresponding to £3 = ^2, which 
should give the radius of the stable circular path, is such that < rs = r2 < r_ = r s i ng . 

5.2.3 One real root for w(p) = 

This case corresponds to (Eq. (13.241) ) 

p4<P<fc (5.26) 

leading to 

Po < e r < p„o < p + < p_ = +00 (5.27) 

where e r is the real root given by (13.25I ). For a < —2 it is not possible to have g2 = and gj = 0, so there is no 
solution e r = with multiplicity 3. In the generic case j3i < j8 < and ^3 7^ there is an absorbed path from 
Poo to the singularity at p = p_ = +00 given by (I3.12ll5.5ll5.8l ). 



6 Conclusion 

To the first order of approximation, all black holes of phantom and normal EMD deflect light in the same 
manner. If we restrict ourselves to physical conditions [a 2 < 1 for T72 = — 1 and a 1 < (1 + y)/(2y) for T72 = + 1], 
then 1) for /larger than some 70, which is likely in (—0.2,0) and depending on the parameters of the black hole, 
black holes of E-anti-M-(anti)-D (regardless the sign of 17 1) theory cause light rays to deflect with larger angles 
than black holes of EM-(anti)-D do. The difference in the angles and the relative discrepancy ever increase with 
a 2 for fixed (u n , 7). For fixed (a 2 , 7), they increase with l/r„ and diverge as r n approaches the photon sphere of 
E-anti-M-(anti)-D black holes. 2) For — 1 < 7 < 70 light is more deflected by the black holes of EMD than by 
those of E-anti-MD, the relative discrepancy for larger values of the impact parameter is, however, much larger 
for the black holes of E-anti-MD. 
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Relativistic images IT331 is another ingredient, besides deflection, allowing for the determination of the 
nature of matter. These observations happen on the photon sphere of the solution. If the latter is different from 
a Schwarzschild black hole, the images are observed for a radial distance larger or smaller than 3M. 

The method based on Weierstrass polynomial to integrate geodesic motion and determine exact solutions is 
equivalent to other methods using potential barriers and can be applied systematically. The advantage of using 
the method based on Weierstrass polynomial is that motion is allowed in at most two regions: in between the 
smallest root of the polynomial and the intermediate one and/or for values gritter than the largest root. This 
highly simplifies the problem. Some of the systematic resolutions consist in: 1) The angle of deflection has a 
standard formula for all problems that can be brought to Weierstrass differential equation. 2) If the two smallest 
and intermediate roots of the Weierstrass polynomial are equal for some value of the vector of parameters, there 
should be a stable circular path (photon sphere) for the corresponding radial coordinate r if the latter is within 
accessible limits to observers. 3) If the two largest and intermediate roots are equal for some value of the vector 
of parameters, there should be an unstable circular path for the corresponding radial coordinate r if the latter is 
within accessible limits to observers. 4) Existence of spiral paths, which approach endlessly the photon spheres, 
is a consequence of any of 1) or 2) or both. 5) Existence and identification of divergencies for the angle of 
deflection: a logarithmic one if 3) holds or a power law one (to the power —1/2) if the three real roots are zero. 
6) Ordering of the parameters expressing spatial infinity, singularity, horizons and so on on the Weierstrass 
axis is derived by circular rotation (from their given order relations on the r-axis) in the one or the other way 
depending on the coordinate transformation relating the Weierstrass radial coordinate to the spherical radial one 
(increasing or decreasing). 

Phantom Reissner-Nordstrom black holes are characterized by the existence of trapped and absorbed null 
paths which do not exist for normal Reissner-Nordstrom black holes. Their other non-common paths include 
many-world (periodic bound) and two-world paths which exist only for normal Reissner-Nordstrom black holes. 
Their common paths include scattering, spiral (existence of logarithmic divergencies) and unstable circular paths 
with radii approaching the Schwarzschild limit from above for phantom black holes and from below for normal 
ones. 

Both phantom cosh and normal sinh black holes of EMD theory are characterized by the presence of scatter- 
ing, trapped and unstable circular paths, thus spiral paths and existence of logarithmic divergencies. The photon 
spheres are larger or smaller than the Schwarzschild one, respectively, and approach it in the limit of no electric 
charge. The phantom solution has no absorbed path while the normal one does have. 
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(a) :Phanton (b) :Normal 




Figure 9: Plots of Y = t — u 2 + 2Mu 3 — rj2q 2 u 4 . (a): The phantom Reissner-Nordstrom black hole(r72 = —1). u r is the lowest 
root of t — u 2 + 2Mk 3 + q 2 u 4 = with «_ < u r . (b): The normal Reissner-Nordstrom black hole(7]2 = +1). u r is the largest root of 
I — u 2 +2Mm 3 — q 2 u 4 = with m_ < M r . For both plots, (u+, «_) are the intersections of the graphs of F = — u 2 + 2Mu 3 — rj2q 2 u 4 with 
the M-axis, which are the same graphs as those shown here but shifted downward I units. 

Appendix: Order relations for the phantom and normal Reissner-Nordstrom 
black holes 

The phantom case Tfe = — 1 

In the case where all the four roots of £ — u 2 + 2Mu 3 + 2q 2 u 4 = have multiplicity 1 we can choose any root to 
perform the reduction of ( 14.31 ) to dl.2| ). For the phantom Reissner-Nordstrom black hole we choose u r to be the 
smallest root as shown in Figure [9] (a) 

m_ < u~ < u+ < M3 < < Ui_ < U\ < u + < +00. (A.l) 

As defined in the first expressions of Eqs. (14.6114.71 ), (C\,Cq) are proportional to the first and second derivative 
of I — u 2 + 2Mu 3 +q 2 u at u = u r , respectively. At the point M r (u r ,0) of the (a)-plot the function is decreasing 
and concave up (convex), thus C\ < 0, Cj > and p is an increasing function of u (a decreasing function of r). 

Since C\ < 0, the coordinate transformation (14.51 ) "splits" the point u r into m+ and u~ , which correspond to 
p = — 00 and p = +00, respectively. If one starts to move on the w-axis from the right to the left: from u = +°° 
(r = 0) to u + to u\ to • • • to m+ to u~ and finally to m_. Since p is an increasing function of u, the corresponding 
point on the p-axis starts to move from po to p+ to e\ to • • • to p = —00 then back in a circular rotation to 
p = +00 and finally to p . Thus, we have the following order relations for the p -parameters 

- 00 < e 3 < po, < e 2 < e\ < p+ < p < p_ < p r = +°° . (A.2) 

[Of course, this ordering can be derived by algebraic methods]. 

If I = (M2 on the w-axis), then u r = w_, U2 = M3 = and u\ = u + (p r = +°°, £2 = £3 and e\ = p + ). 
Using (T3~T4l 1431 K1\WI[ K9\ 02114031) we obtain 12p = 5 - [6772(1 + y/ 1 - fya 2 ) / a 2 ] (with tj 2 = -1) and 
the latter equation becomes 

ej, = e 2 = Poo = --^r < e\ = p+ = \ < p < p_ = . (A.3) 
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If £ = £ + , then u\ = u 2 (Mi on the w-axis) and (IA.2I) becomes 



- 00 < e 3 < p„ < e 2 = e\ < p+ < p < P- < p,- = +°° . (A.4) 
If £ = then u r = 113 (M3 on the i<-axis) and the root u r has multiplicity 2. 

The normal case t]i = +1 

In the case where all the four roots of I — u 2 + 2Mm 3 — 2q 2 u 4 = have multiplicity 1 we can choose any root to 
perform the reduction of ( 14.31 ) to dl.2| ). For the phantom Reissner-Nordstrom black hole we choose u r to be the 
largest root as shown in Figure [9] (b) 

W3 < < U 2 < U\ < U + < U- < U~ < < +00 . (A.5) 

At the point M r (u r ,0) of the (b)-plot the function is decreasing and concave down (concave), thus C\ < 0, 
C2 < and p is an increasing function of u (a decreasing function of r). Since C\ < 0, the coordinate transfor- 
mation (14.51 ) "splits" the point u r into and m,7, which correspond to p = —00 and p = +00, respectively. If 
one starts to move on the w-axis from the right to the left: from u = +°° (r = 0) to m+ to m,7 to m_ to ■ • • to u 2 
to W3 and finally to u = —00. Since p is an increasing function of u, the corresponding point on the p-axis starts 
to move from po to p = —00 then back in a circular rotation to p = +00 to p to • • • to e 2 to es and finally to po 
again. Thus, we have the following order relations for the p -parameters 

- 00 < p < e 3 < p„o < e 2 < e\ < p + < p_ < p,- = +00 . (A.6) 

If I = (M 2 on the w-axis), then u r = U-, u 2 = u 3 = and u\ = u+ (p r = +00, e 2 = e 3 and e\ = p + ). 
Using (l3^l4310l4T^l4^l4Tl^l4TT3l) we obtain 12p = 5 - [6tj 2 (1 + V 1 = ma 2 ) /a 2 ] (with tj 2 = +1) and 
the latter equation becomes 

- 00 < p < e 3 = e 2 = poo = -^r < ei = p+ = ^ < P- = +°° (A.7) 
If £ = then wi = «2 (A/3 on the w-axis) and (IA.6I ) becomes 

- 00 < po < e 3 < poo < e 2 = e x < p + < p_ < p r = +00 . (A.8) 

If ^ = < 0, there is still one real root u + < u r = u\ < «_ (Mi on the w-axis). This case is excluded, 
however, if it were possible for a photon to move with a negative energy, it could do it on a confined stable circle 
with radius r_ < r = M(3 — V9 — 8a 2 ) /2 < r + which shrinks to zero as a 2 = q 2 /M 2 approaches zero. 
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